Spin precession and alternating spin polarization in spin-3/2 hole systems 
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The spin density matrix for spin-3/2 hole systems can be decomposed into a sequence of multi- 
poles which has important higher-order contributions beyond the ones known for electron systems 
[R. Winkler, Phys. Rev. B 70, 125301 (2004)]. We show here that the hole spin polarization and 
the higher-order multipoles can precess due to the spin-orbit coupling in the valence band, yet in 
the absence of external or effective magnetic fields. Hole spin precession is important in the context 
of spin relaxation and offers the possibility of new device applications. We discuss this precession 
in the context of recent experiments and suggest a related experimental setup in which hole spin 
precession gives rise to an alternating spin polarization. 

PACS numbers: 72.25.Dc, 71.70.Ej, 72.10.-d, 73.23.-b 



Spin electronics is a quickly developing research area 
that has yielded considerable new physics and the 
promise of novel applications Q. Among the main 
focuses of spin electronics are semiconductor systems, 
where ferromagnetic semiconductors and spin polarized 
transport stand out as major areas of interest. In these 
fields, the importance of holes as compared with elec- 
trons is manifold. Firstly, the compound semiconductors 
exhibiting ferromagnetism are p-type materials, in which 
ferromagnetism has been shown to be mediated by the 
itinerant holes 01 ■ These materials have also been used 
as sources of spin-polarized holes, motivating the search 
for a more complete understanding of hole spin preces- 
sion and relaxation. Secondly, the spin-Hall effect was 
first studied in the context of a hole Hamiltonian 0, 0| 
and first observed in a two-dimensional hole gas . With 
these facts in mind, we devote this Letter to an in-depth 
study of the spin dynamics of hole systems. 

The conduction band of common semiconductors like 
GaAs is described by a spin-1/2 Hamiltonian so that elec- 
tron spin dynamics are relatively amenable to theoretical 
investigation. It is well known that for spin-1/2 electrons 
the spin-orbit interaction can always be written as a wave 
vector-dependent effective magnetic field 0, . Electron 
spin precession in this effective field, as well as in an ex- 
ternal field, has been discussed extensively and is well 
understood The spin dynamics of holes, described 
by an effective spin s = 3/2 0, have been studied to a 
lesser extent [lOllllI- In general, the coupling of the hole 
spin and orbital degrees of freedom cannot be written as 
an effective magnetic field, so the simple picture of a spin 
precessing around an effective Zeeman field is not appli- 
cable to holes. Nevertheless, we will show that the spin 
dynamics of hole systems can be viewed as a precession, 
if precession is understood as a nontrivial periodic mo- 
tion in spin space described by an equation of the type 
^ = ji[Tl,S] for a suitably generahzed spin operator S 
and spin Hamiltonian Ti. While it has long been known 



that holes lose their spin information much faster than 
electrons, to our knowledge a quantitative discussion of 
hole spin precession based on a firm theoretical footing 
has not been attempted in the literature. In particular, 
we will show that in appropriate geometries the spin po- 
larization can alternate in time, even at magnetic field 
B = 0. This alternating spin polarization is a dynamical 
effect specific to spin-3/2 systems. 

The explicit form of the spin-orbit interaction depends 
sensitively upon the symmetry of the system |l2l 1131 |lj| , 
a fact which is captured in an invariant decomposition 
of the spin density matrix [Tsf . Neglecting small terms 
with cubic symmetry, the spin density matrix of spin-3/2 
systems can be expressed in terms of invariants of SU(2). 
These are proportional to spherical tensors - a monopole, 
dipole, quadrupole and octupole, denoted respectively by 
A/" (a multiple of the identity matrix), (the vector 
of spin- 3/2 matrices), and M^: 



P - 



(1) 



The dot product for spherical tensors appearing in Eq. 

is defined in Ref. [ig. The moments po, <S, Q, and 
O provide a set of independent parameters characterizing 
the weights of the multipoles in the expansion [isj . 
Excepting normalization factors, the monopole moment 
Pa is the carrier density, while the dipole moment <S corre- 
sponds to the spin polarization (Bloch vector) at B > 0. 
The quadrupole moment Q reflects the splitting between 
the heavy holes (HHs, spin-z projection rUg = ±3/2) and 
the light holes (LHs, = ±1/2). The octupole moment 
O is a unique feature of s = 3/2 systems at i3 > 0. We 
note that in the density matrix of spin-1/2 electrons only 
the first two terms of the expansion ^ are present. We 
would like to point out that the term "multipole" is used 
here in a rather different sense from that of Ref. ^ . 

The dynamics of spin-3/2 hole systems are determined 
by the 4x4 Luttinger Hamiltonian H |^, , which, in 
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the spherical approximation, reads 
^'71 ^0 
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where 71, 7, and /c are Luttinger parameters, Too is the 
bare electron mass, and are tensor operators fisf . 
In the lowest order of k and B we have /C" = fc^ and 
K} = B. The operator K,^ is responsible for the HH-LH 
coupling in hole systems. We neglect the terms with cu- 
bic symmetry, which are small corrections in the present 
context, as discussed in Refs. [isllr^ . 

We will analyze in detail the Heisenberg equations of 
motion (HEM) of the multipoles Equation lO sug- 
gests that we study first the HEM 
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The HEM of a dipole TVi'^ evolving under the action of the 
Zeeman term Ti} describe the well-known Larmor preces- 
sion of holes in a magnetic field . These equations are 
closed, i.e., the RHS does not depend on the other mo- 
ments , j" 7^ 1. In general, however, the HEM for 
the in spin-3/2 hole systems cannot be decoupled, 
having important consequences for hole spin precession. 
This can be seen from Tabled which gives the invariant 
decomposition of the RHS of Eq. ((2Jl for different j and 
j' The most remarkable entry in the table is the 

one for the quadrupole iW^ propagating in time due to 
a quadrupole term in the Hamiltonian, H^, which rep- 
resents the spin-orbit interaction that gives rise to the 
HH-LH coupling. The HEM for are not closed. A 
quadrupole M"^ precessing in a quadrupole field "decays" 
into a dipole and an octupole . This implies 
that spin precession of an initially unpolarized system 
can give rise to spin polarization, even though there is no 
external or effective magnetic field. [As mentioned above, 
we use the term spin precession for any HEM Q with 
> 0.] Moreover, spin precession of hole systems does 
not preserve the magnitude of the Bloch vector, as will 
be shown below explicitly. These observations point to 
the fact that the spin dynamics of spin-3/2 hole systems 
are qualitatively different from those of spin-1/2 electron 
systems, where the HEM for the monopole (number den- 
sity) is decoupled from that of the dipole (Bloch vector 
(t). As a result, electron spin precession preserves the 
length of the Bloch vector (cr), i.e., d\{(T)\/dt ~ 0. 

We want to determine and interpret the explicit time 
evolution of p in order to investigate the physics ab- 
sent in electron systems. This calculation is most eas- 
ily carried out in the Schrodinger picture, which reflects 
the equivalence of the Heisenberg and Schrodinger pic- 
tures for this problem. In the absence of external fields 
and disorder, the density matrix satisfies the quantum 
Liouville equation ^ = -I [p, 7i] , with formal solution 



TABLE I: Irreducible representations T>-' of SU{2) of the 
(linear combinations of) multipoles M'' contained in an in- 
variant decomposition of [l/ih) [M-' ]. 
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p{t) = e"''"*/'^p(0)e''"*/'^, where e*^*/'* is the time evo- 
lution operator (which can often be evaluated in closed 
form) . We consider first an example where we assume the 
hole spin to be oriented initially along the z-direction, 
TOs = +3/2, so that 
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where the tilde indicates that p has been normalized with 
respect to the total density 2po- This equation demon- 
strates that, in general, the density matrix of holes can- 
not be written simply as the sum of a monopole and a 
dipole. The higher multipoles will be present We 
want to restrict ourselves to H = TiP + 0.^, i.e., B = 0. 
Table|l]shows that p evolves as a combination of a dipole, 
a quadrupole, and an octupole. The implications of this 
fact for the spin are seen by following the motion of the 



Bloch vector S{t) 
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S{t) 



Sit) -= ^{ So [cos2(tji) +c4sin2(tjt)] 

+ fees (1 + c^) sin^([ji) (5) 
-|- iSo X fc) 2 c s sin(ijji) cos(a;t) } , 

where the unit vector Sq is the orientation of S' at t = 0; 
u! = '^hk'^/mQ = (e/j — ei)/2h, with e/j and e; the HH 
and LH energies, respectively; c = S'o ■ fc is the cosine 
of the angle between So and fc; and s is the sine of the 
same angle. When So is parallel to fc (i.e., c = 1), we get 
Sit) ~ So, which is due to the fact that the initial state 
is an eigenstate of the Hamiltonian. In general, neither 
the magnitude nor the orientation of the Bloch vector 
are conserved. This is illustrated in Fig. ^ showing Sit) 
for an angle of 60° between Sq and fc. Helicity S ■ k 
is conserved, a well-known fact about this model, which 
sheds additional light on spin precession in hole systems. 
Since -^iS • k) = and the wave vector is not changing, 
^•fc = 0. Therefore, whenever the magnitude of the spin 
changes, the angle between spin and wave vector must 
change in order to preserve the projection of S onto fc. 
As a consequence, no nontrivial spin precession occurs 
when 5 _L fc (or 5' II fc). We note that energy is also 
conserved for spin precession in hole systems. Yet for 
B > 0, when H = n^ + + Ti'^, it can be shown that. 




FIG. 1: Precession of the Bloch vector S (bold arrow) about 
the wave vector k when the angle between the two is initially 
60°. We have assumed S{t = 0) is parallel to the z axis, and 
k lies in the xz plane. 



in the general case, energy is transferred back and forth 
between and Ti^ as time progresses. 

We would like to relate our work to recent experi- 
ments by Grayson et al. [20I [2l| which demonstrated 
that a high-quality bent heterostructure can be grown 
on top of a pre-cleaved corner substrate that allows one 
to drive the charge carriers around an atomically sharp 
90° corner. Grayson's experiments were performed on a 
two-dimensional (2D) electron system in a GaAs/AlGaAs 
quantum well. We show here that a similar system con- 
taining holes gives rise to fascinating new physics. The 
setup is sketched in Fig. [2Ia) . We assume B = so that 
the Hamiltonian is H = -I- H^. An unpolarized HH 
wave packet travels in the 2D channel Li in the +x di- 
rection. For simplicity we assume that all orbital degrees 
of freedom are integrated over, so that p depends only on 
the spin indices and does not depend on the wave vector 
fc. Nevertheless, the coupling of spin and k is present 
in the Hamiltonian through Ti.^. The magnitudes of the 
normalized moments at i < (i.e., before reaching the 
corner) are S ^ O = Q and Q = 1/2. The spin quan- 
tization axis of the HH states in Li is parallel to the 
0-direction. After the wave packet has passed the cor- 
ner, the HH states are not eigenstates of Ti,^, their spin 
quantization axis being perpendicular to the spin quan- 
tization axis supported by the quadrupole field in L2. 
Therefore, the quadrupole and octupole moments in L2 
oscillate in time 

|Q(i)|2 = (1/16) -h (3/16) cos2(tj,t) (6a) 
|d(i)|2 = (3/16)sin2(a;,t), (6b) 

with precession frequency Uz — {£h — £0/2^ — 
/ rriQuP' . This frequency can be tailored by vary- 
ing the width w of the 2D channel. For GaAs in the 
spherical approximation we have 7 = 2.58 and we take 
w = IQ nm, yielding a precession period 2t: /lOz ~ 11 ps. 

For the simplified geometry of Fig. Ela) the spin po- 
larization S in L2 remains zero, as required by the con- 



FIG. 2: A bent structure allows holes to be driven around a 
corner [2(ll2H. (a) An idealized representation of the corner 
device, (b) A more realistic model of the corner, as discussed 
in the text. In regions C and L2 the spin polarization oscillates 
as a function of time. 



servation of helicity discussed above. However, assuming 
a sharp 90° corner for an ideal 2D system is certainly an 
oversimplification, even when the corresponding quasi 2D 
system has atomically sharp interfaces 20] . A more real- 
istic treatment can be obtained by modeling the transi- 
tion region between the channels Li and L2 as a sequence 
of two 45° corners as sketched in Fig.l^b). Once again, 
an unpolarized HH wave packet travels in channel Li in 
the +x direction, the initial conditions being the same as 
in the previous example. If the HH wave packet enters 
the central region C at t = 0, we get for the squared 
normalized moments in this region 

|5(i)p = (9/80) sin^{ujzt) (7a) 
|Q(i)|2 = 1/64 (15/64) cos2(a;^t) (7b) 
\0{t)\^ = (39/320) sm^{ujzt). (7c) 

Equation IjTall shows that the initially unpolarized hole 
current acquires an alternating spin polarization S{t) 
due to spin precession at S = 0. In Cartesian co- 
ordinates, the Bloch vector in region C reads S{t) = 
[0, — s\n{ijJzt),Q\. When the HH wave packet enters 
the channel L2 it continues to precess. We get for the 
spin polarization 
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Here T is the time required to traverse C which depends 
on the length of C and the magnitude of the in-plane wave 
vector. It determines the fraction of the spin polarization 
in L2 that will be oscillating. If we take the length of C 
to be of the order of the channel width, namely w = 
10 nm, and the initial wave vector kp = 0.1 nm~^, we 
get T ~ 0.5 ps. The amplitude of Sy in this case is 
approximately 0.1. We omit here the qualitatively similar 
but more complicated expressions for Q,{t) and Oit). We 
note also that the approach in Fig. [2Ib) can be further 
extended in a transfer-matrix-like approach in order to 
describe more complicated geometries. 

We expect hole spin precession to be robust against dis- 
order provided lOz > I/t's, a range experiment suggests 
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is realistic. Hole spin relaxation times r, in the range 
of tens of ps ^ or even one ns have been observed 
experimentally for GaAs quantum wells, demonstrating 
that the spin precession discussed in this paper cannot 
be neglected in the context of hole spin dynamics. We 
expect that the precession can be observed experimen- 
tally using time-resolved optical techniques such as Kerr 
or Faraday rotation spectroscopy. 

The examples above identify an interesting potential 
application of spin-orbit effects in hole systems. By 
choosing appropriate initial conditions the spin polariza- 
tion in the bulk of the system can be made to oscillate as 
a function of time. This oscillating spin polarization is to 
be contrasted with the Zitterbewegung discussed lately 
j24l | , in which the carrier position oscillates as a function 
of time. Unlike conventional transistors, which work in 
the diffusive limit, the examples we discuss are valid in 
the ballistic limit, a feature shared with the Datta spin 
transistor . The situation we describe is nevertheless 
different from the spin precession in a Datta spin transis- 
tor where the electrons precess not as a function of time 
but as a function of position |^ . 

The above example illustrates the controlled hole spin 
precession in a confined geometry in a ballistic regime, 
where the holes experience an essentially unidirectional 
quadrupole field. The situation is qualitatively different 
in, e.g., bulk material in a diffusive regime, where the 
occupied hole states are characterized by many different 
orientations of the wave vector k. The quadrupole field 
is thus randomly oriented so that spin precession in this 
field must also be taken into account for spin relaxation, 
similar to Dyakonov-Perel (DP) spin relaxation in elec- 
tron systems. Previously |2(l l27l|. spin relaxation of hole 
systems was discussed mostly using an Elliot- Yafet-type 
approach , where spin relaxation is caused by momen- 
tum scattering events. DP spin relaxation, on the other 
hand, implies that the spin orientation is lost due to spin 
precession in between the momentum scattering events. 
Such spin precession for the case of holes has been de- 
scribed by Averkiev et al. but only as a result of 
the spin splitting induced by bulk and structure inver- 
sion asymmetry in 2D systems. On the other hand, the 
spin precession described by Eq. Q is different in na- 
ture, being due to the quadrupole field so that it is 
present also in inversion symmetric systems. Note also 
that hole spin precession due to is typically one or sev- 
eral orders of magnitude faster than the precession in the 
effective magnetic field due to a broken inversion symme- 
try so that Ti.^ can yield the more important contribution 
to DP hole spin relaxation. (We use here the term DP 
spin relaxation in a generalized sense for any relaxation 
based on spin precession, even if motional narrowing 7| 
is not important.) Equation © suggests that hole DP 
spin relaxation in bulk systems occurs on time scales of 
the order of 2nu)~^ — 2T:{2^hk^ /mo)~^ . In bulk GaAs, 
for holes optically excited with light at 800 nm, recent 



experiments yielded spin relaxation times of about 110 fs 
|29l |. while the simple estimate above gives ^ 200 fs for 
this case. (This number depends sensitively on the wave 
vector of the occupied hole states that can only be esti- 
mated for the experiments in Ref. [1^.) We remark that, 
in analo gy with spin relaxation in anisotropic electron 
systems |28t |30| , the spin relaxation of hole systems will 
be characterized by different relaxation times for each 
multipole moment in the density matrix. 
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